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Abstract. We argue that the well known Dynes formula [Dynes K C et al. 1978 
Phys. Rev. Lett. 41 1509] for the superconducting quasiparticle density of states, 
which tries to incorporate the lifetime broadening in an approximate way, cannot 
be justified microscopically for conventional superconductors. Instead, we propose a 
new simple formula in which the energy gap has a finite imaginary part — A2 and the 
quasiparticle energy is real. We prove that in the quasiparticle approximation 2A2 gives 
the quasiparticle decay rate at the gap edge for conventional superconductors. This 
conclusion does not depend on the nature of interactions that cause the quasiparticle 
decay. The new formula is tested on the case of a strong coupling superconductor 
Pbo.gBio.i and an excellent agreement with theoretical predictions is obtained. While 
both the Dynes formula and the one proposed in this work give good fits and fit 
parameters for Pbo.gBio.i, only the latter formula can be justified microscopically. 
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Almost thirty years ago Dynes, Narayanamurti and Garno [Ij proposed that the 
quasiparticle recombination time in a strong-coupled superconductor can be directly 
measured from the width of the peak in the tunneling conductance dI{y)/dV of a 
superconductor-insulator-superconductor tunnel junction at the sum of the gaps. They 
found that the data on Pbo.gBio.i-insulator-Pbo.gBio.i planar tunnel junction could be 
fitted quite well for voltages near twice the gap if the quasiparticle density of states 

p{E) = Re , ^ (1) 

Ie^-/s?{e) 



in the expression for the tunneling current 

/ + 00 
dEp{E)p{E + eV)[f{E) - f{E + eV)] (2) 
-oo 



is replaced by 



PoiE, Yd) = Re , ^ '^"^ , (3) 
'{E - zr^)2 - Al 



with real and E-independent and the measured gap edge Aq. In ([T]) A{E) is the 
complex gap function and / and e in ([2]) are the Fermi function at temperature T and 
the magnitude of electron charge, respectively. It was proposed P that the temperature 
dependent parameter in ([3]) incorporates the quasiparticle lifetime effects. A good 
agreement between the measured F/)(T) and a microscopic calculation [1] based on 
the work by Kaplan et al. [2] for a number of temperatures below the transition 
temperature of Pbo.gBio.i was taken as a justification for the replacement of p{E) with 
Pd{E, Yd) and for the interpretation of parameter 2Fd as the inverse of the quasiparticle 
recombination lifetime. Formula ([3]) is now widely known as the Dynes formula and it 
has been applied to a variety of low temperature (T <C Tc) tunneling experiments 
ranging from tunneling into the bulk [3] and thin film P| inhomogeneous/granular 
superconductors to the tunneling into a two-band superconductor MgB2 [5] and 
tunneling into a novel superconductor CaCg P, [7]. The Dynes formula was also 
recently used to describe the density of states obtained in photoemission studies of 
superconducting h-ZrRuP [8] and of filled skutterudite superconductor LaRu4Pi2 [9]. 

However, the ansatz cannot be justified for a conventional strong coupling 
superconductor, such as Pbo.gBio.i [I], from first principles. Indeed, p{E) is given in 
terms of the diagonal electron Green's function in the superconducting state 

^> - E^Z^ik,E)-<P^ik,E)-el ' 
where Z is the complex renormalization function and is the complex pairing self-energy 
[iniin], as 

p(E) = -^^^lra^Gn{KE), (5) 

where N{0) is the normal state density of states at the Fermi level. All interactions enter 
via the self-energy terms Z and 4> and assuming that they do not depend on momentum 
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p{E) = Re ^ ^^(^^ (6) 



k one finds 



= Re , ^ (7) 
^E^-A^{E) 

where in the last step Z{E) and (f){E) have been eliminated in favor of the gap function 
A[E) = (j)[E)/Z[E). Clearly, all the lifetime effects which enter via (piE) and Z{E) are 
ultimately incorporated in the complex gap function A{E) and the tunneling current 
I{V) depends on the full complex gap function as is clear from equations ([1]) and ([2]). 
Note that ([6]) cannot be cast into the form ([3]) by a suitable choice of Z{E) (e.g. taking 
Z{E) = 1 - iTn/E would give Re [{E - iTD)/^{E - iVoY - where the pairing 

self-energy appears instead of the gap A, and the measured dl/dV gives A and not 
0). 

Instead of replacing p{E) with Pd{E, Yd) it is more reasonable to keep A{E) in ([T]) 
constant but complex for E not too far from the gap edge Aq, i.e. replace ([1]) with 

PA A2) = Re , ^ , (8) 

_ (Ao - ZA2)2 



where — A2 is the imaginary part of the gap aX E = Aq. It is well known that at a 
finite temperature the imaginary part of the gap at the gap edge is finite as a result 
of quasiparticle damping (see figure 45 in [H]). In fact, it is easy to prove that in the 
quasiparticle approximation [2] the quasiparticle decay rate at the gap edge is equal to 
-2ImA(ii^ = Aq). Assuming that at ii^ = Aq the imaginary parts Z2 and and 02 of Z 
and 0, respectively, are much smaller than the corresponding real parts one finds that 

-ImA(E = A„)«^2^^1^MzM£^, (9) 

^ Zi(0) ^ ' 

where Zi{0) is the real part of Z{E = 0). Expression ([9]) is identical to the equation of 
Kaplan et al. for the quasiparticle decay rate parameter V{E = Aq) [2] (see equation 
(5) in [2J). This result is quite general and does not depend on the specific interactions 
leading to quasiparticle damping, i.e. whether it is the electron-phonon interaction 
which was considered in [H [2] , or the dynamically screened Coulomb interaction in the 
presence of disorder which was assumed to be the cause of lifetime broadening in low 
temperature tunneling experiments into three-dimensional granular aluminum [3] and 
quench-condensed two-dimensional films of Pb and Sn [4t\ . All that is required for 

2r(k, E = Ao) = -2 Im A(k, E = Aq) , (10) 

to be valid, where 2r(k, = Aq) is the inverse quasiparticle lifetime with k on the 
Fermi surface, is that the imaginary parts of 0(k, E) and Z{k, E) are much smaller 
than their respective real parts near the gap-edge. Needless to say, ffTOj) does not apply 
to unconventional superconductors characterized by X)kGF5^(k) =0, where FS is the 
Fermi surface, for k near the gap nodes |12j . 
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In the case of Pbo.gBio.i we find that equation ^ produces fits to dI{y)ldV which 
are at least as good as those obtained with the Dynes formula (j3]). Instead of trying to 
fit the original data from [1], which in addition to the temperature dependent lifetime 
broadening were assumed to contain an intrinsic (background) width of O.OlmeV, we 
fitted dI{y)/dV calculated from the solutions A(i?) and Z{E) of the finite temperature 
Eliashberg equations [IHl [H] on the real axis using the Eliashberg function 
for Pbo.gBio.i [I3]. Thus, the width of the peak in our calculated dI{y)/dV arises solely 
from the temperature dependent lifetime broadening and we could compare directly the 
value of the fit parameter A2 in equation ([8]) with our solution — \m.A{E) for E at the 
gap edge. Moreover, we could calculate the decay rate parameter T{E) directly from 
our solutions of Eliashberg equations [2j (see equation (4) in |2]) 

T{E) = EZ,{E)/Z,{E) - ME)ME)/[Z',{E)E] (11) 

and compare its value aX E = Aq with A2 obtained from the fits with equation ([8]). 
We note, however, that there is a good agreement between the shapes of the calculated 
dI{V)/dV and the measured ones p^J down to T =2.75 K as illustrated in figure [1] 
for T=3.5 K. In figure [1] the results are plotted as functions of eV — 2Aq since with 
our choice of the Coulomb pseudopotential /i*(a;c)=0.1034, which was fitted to the 
experimental zero temperature gap edge Ao=1.54 meV [13| for the cutoff a;c=100 meV 
in the Eliashberg equations, we obtain somewhat higher values of Ao than those found 
in [1]. As the Coulomb pseudopotential term in the Eliashberg equations is purely real 
it does not affect the imaginary parts of the solutions [lOl [11] . 

In figure [2] we show the fits to the calculated dl/dV using the Dynes formula 
([3]) and the formula with the complex gap (IHl). On the scale of figure Ej which was 
chosen to match the scale of figure 2 in [Ij, both equations ([3]) and ([8]) give equally 
good fits. Moreover, the values of the fit parameter A2 turn out to be nearly the 
same as the values of the fit parameter F^) at all temperatures considered. One can 
understand why two different functional forms and ([H]) give nearly identical fits to 
dl/dV with nearly identical fit parameters ^ A2 from the fact that in the limit 
Yd, A2 <C Ao the approximations ([3]) and ([8]) to p(Ao) give Aq/Td/2 and AqI A2/2, 
respectively and the height of the peak in dl/dV is most sensitive to the maximum in 
the quasiparticle density of states. However, it is clear that as the lifetime broadening 
grows compared to the gap edge the difference between the fit parameters obtained with 
([3]) and with ([8]) increases and the quality of fits with the Dynes formula deteriorates 
compared to the fits with ([8]) as illustrated in figure [3], in particular at lower voltages. 
The reason is that for F^i, A2 -C Ao in the limit of small energy pd{E,To) = To/Aq, 
while Pj\{E, A2) = {A2/Al)E to the first order in E, i.e. p£){E,T£)) does not vanish 
at E=0. We note that the experimental low-temperature densities of states obtained 
for three-dimensional granular aluminum [3] do vanish at £"=0 (see figure 3 in [3J) , 
while those obtained for two-dimensional quench- condensed tin films [4j do not (see 
figure 2 in [1]). The precise reason for such a difference between three-dimensional and 
two-dimensional disordered conventional superconductors is not known at the present 
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time. 

As one could have expected, the fitted values of A2 turned out to be equal to 
the imaginary parts of our solutions A(i?) of the Eliashberg equations at -E = Aq to 
within a few percent at all temperatures considered. The values of A2 extracted from 
the fits to the calculated dl/dV agree with the values of the fit parameter r(= Yd) 
reported in [1] before correction for the background to within a percent or two down to 
T=4.2 K. At T=3.5 K the difference is about 30% and yet the shapes of the calculated 
and measured dl / dV in figure [1] seem to agree quite well. A further reduction of the 
measured F by the background value of O.OlmeV would increase the difference between 
the lifetime broadening parameters to about 150%. At T=2.75 K our fitted value (the 
fit is not shown here) is A2=0. 00226 meV which is 80% lower than the measured F [T] or 
more than twice the measured value after the correction for the background. It is quite 
plausible that at low temperatures, when both the experimental and the theoretical 
data in the peak change very rapidly, it is difficult to determine the actual maximum 
in dl/dV to which the fit parameters are most sensitive. It is likely that the maximum 
in dl/dV gets underestimated at low T having as a consequence too high values of 
the lifetime broadening parameter. We believe that is the reason for the discrepancies 
between our fitted values of A2 and those found in [I] at low temperatures and that 
there is no need to invoke the intrinsic temperature-independent broadening parameter. 

Finally, in figure H] we show the temperature dependence of the quasiparticle lifetime 
T at the gap edge obtained from h/r = 2A2 (open squares) and ^/r=2(F-0.01meV) 
(filled circles) with the values of F taken from figure 2 in . In the same figure we show 
theoretical predictions for the recombination time (solid line) and the total lifetime 
r (dashed line) at the gap edge based on approximate equations of Kaplan et al. [2] 



n 

— = c dna\n)F{n) 

Tr J2A0 



2Ao ^(Q - Ao)2 - Al 

x^^^^Jn{Q) + l]fiQ-Ao), (12) 

^ ^ ,^ 9.^^ n + A 



- = c rdna\n)F{n) , "^^q ^ 

^« -^0 . /(fi + Ao)2 - Al 



•■^_^^^M^)[l-/(^ + Ao)], (13) 

where n{Q) is the Bose function, C = 27r/{Zi(0)[l — /(Aq)]} and h/r = h/rr + h/Ts. 
A good agreement between the measured ^/2(F-0.01meV) and calculated according 
to (IT^ was taken as a justification of the Dynes formula in |lj. We note that the 
integrand in (fT2|) has a square root singularity at the lower limit of integration which 
has to be handled analytically if Tr is not to be overestimated. Comparing figure 3 in [1] 
and figure 4 in this work it is clear that our calculated Tr is considerably lower at the low 
temperatures than the one calculated in [1] as the filled circles in both figures represent 
/i/2(F-0.01meV). In addition, we show in figurelHthe lifetime calculated directly from our 
solutions of the Eliashberg equations in the quasiparticle approximation H/t = 2F(Ao) 
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(plus signs), where r(Ao) is computed using equation (ITT]) . The agreement between the 
values for the total quasiparticle lifetime r at the gap edge obtained from the fits with 
formula (IHD and both theoretical predictions is excellent. 

In conclusion, we have shown that one can, indeed, obtain the total quasiparticle 
lifetime at the gap edge from the fits of the derivatives of the I — V characteristic of 
a superconductor- insulator-superconductor tunnel junctions using equation ([8]). The 
interpretation of the parameter 2A2 as the quasiparticle decay rate at the gap edge is 
microscopically justified. While the Dynes formula ([3]) gives correct values for the total 
quasiparticle lifetime, it cannot be justified for conventional superconductors. Hence 
the fact that it works, at least for the cases when the quasiparticle decay rate is less 
than about 20% of the gap edge, is a pure accident. It is likely that for larger values of 
2r/Ao, which seems to be the case in LaRu4Pi2 (2r/Ao ~50%) [9j, equations ([3]) and 
(IE]) would give qualitatively and quantitatively different results. 
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Figure 1. The calculated dl/dV (solid line) and the experimental data points (filled 
circles) from 1, at T—3.5 K. The data are plotted as a fimction of eV-2AQ. 



/O-n -the nriT'Y'or'i' fn/ynmiiln frtr -the li-fo-tnTn c hmn rl en cri cnincr'r'nfn rIomi'O'n n en coi-oi r\-f Q-tn-toc Si 



3.2 



V(mV) 

2.8 2.6 2.4 2.2 



T=3.5K 




A =0.0102meV 



rp=0.0102meV 



A2=0.0254meV 
rp=0.0254meV 



T=5K 



A2=0.0517meV - 
rj^=0.0517meV - 




2.4 2.2 



1.8 1 



T=6.5K A2=0.088meV rj^=0.090meV- 



> 

-a 
-a 



2.2 



1.8 1.6 1.4 1.2 
V(mV) 



Figure 2. The calculated (filled circles) dl/dV at six different temperatures as a 
function of voltage and their fits with ([8]) (solid line) and with the Dynes formula Q 
(dashed line) with A2 and F^i as the only fit parameters, respectively. 
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Figure 3. The calculated dl/dV (dots) at T=6 K versus voltage and the fits with 
JS]) (sohd hue) and ([3]) (dashed line), with A2 and F/) as the only fit parameters, 
respectively. 
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Figure 4. Quasiparticlc lifetime at the gap edge determined in various ways (see the 
text) as a function of inverse temperature. 



